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AN  ANALYSIS  OF  SMALL  SCALE  GAP  TEST  SENSITIVITY  DATA  USING  POROSITY  THEORY  AND 
NONREACTIVE  SHOCK  HUGONIOTS 

The  work  leading  to  this  report  was  done  under  the  task  "Explosion  Initiation 
and  Safety"  SF33-35^-3lU/l8L62.  This  report  permits  the  use  of  small  amounts 
of  experimental  data  at  one  set  of  conditions,  to  generate  critical  sensitivity 
data  over  a broad  spectrum  of  conditions.  This  information  is  of  interest  in 
safety  and  reliability  studies  and  design  of  explosive  systems. 

The  author  acknowledges  the  work  of  Dr.  D.  J.  Pastine  and  Dr.  R.  R.  Bernecker 
in  the  development  of  the  porosity  theory  without  which  the  work  reported  here 
could  not  have  been  done. 
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I.  INTRODUCTION  AND  BACKGROUND 

The  Small  Scale  Gap  Test  (SSGT)1  is  a method  for  indirectly  determining  the 
shock  sen' itivity  of  an  explosive.  In  this  test  a shock  is  transmitted  into  an 
explosive  through  an  attenuating  barrier,  a columr  of  polymethyl -methacrylate 
(PMMA).  The  shock  sensitivity  is  defined  as  that  shock  pressure  required  to 
initiate  the  detonation  of  the  explosive  50$  of  the  time. 

The  SSGT  supplies  information  on  the  shock  strength  needed  in  the  PMMA,  at 
the  PMMA/explosive  interface,  in  order  to  letonate  the  explosive  which  is  at 
some  given  initial  porosity. 

It  would  be  more  meaningful,  however,  to  know,  or  to  be  able  to  calculate, 
the  pressure  in  the  explosive  rather  than  in  the  attenuator.  This  can  be  done 
by  the  usual  impedance  matching  technique  if  the  appropriate  Hugoniot  equation 
of  state  is  known  for  the  explosive  as  well  as  for  the  attenuator.  To  obtain 
a nonreactive  Hugoniot  for  the  explosive,  even  at  a single  porosity,  requires 
considerable  time  and  experimental  effort;  to  measure  Hugoniots  at  a number 
of  porosities  is  ordinarily  economically  unfeasible.  It  is  therefore  desirable 
that  some  method  be  available  to  generate  from  a single  available  nonreactive 
Hugoniot  determined  for  an  explosive  at  a given  density  (porosity),  the  Hugoniots 
at  any  other  stated  densities  for  that  explosive.  It  would  also  be  desirable, 
given  the  sensitivity  of'  an  explosive  at  some  density,  to  be  capable  of  predicting 
the  sensitivity  at  any  othei  stated  density.  It  is  Just  such  a set  of  relationships 
that  this  report  evaluates. 

The  theory,  which  is  detailed  elsewhere,2  provides  a method  for  computing 
the  shock  Hugoniots  of  explosives  as  a function  of  their  initial  porosities, 
provided  one  Hugoniot  at  a stated  porosity  is  known.  This  allows  one  to  compute 
from  SSGT  data,  by  impedance  matching,  the  pressure  transmitted  into  the 
explosive  by  the  PMMA,  since  the  Hugoniot  for  PMMA  is  available.  The 
procedural  technique  for  making  the  calculation  is  detailed  in  Appendix  A.  In 
addition,  the  concept  provides  a set  of  equations  relating  the  shock  pressure 
within  the  explosive  to  the  thermal  energy  density,  E , , immediately  behind  the 

v 

shock  wave.  For  a given  explosive  and  test  conf? duration,  E by  the  concept  is 

X* 

constant  at  the  50$  initiation  point.  The  equations  which  relate  to  tb^ 

evaluation  of  E are  given  in  Apperdix  B.  In  Table  1 are  listed  the  values  of 
_ T* 

*J.  N.  Ayres,  L.  J.  Montesi,  and  R.  J.  Bauer,  "Small  Scale  Gap  Test  (SSGT)  Data 
Compilation:  1959-1972  Volume  I Unclassified  Explosives,"  NOLTR  73-132, 

26  Oct  1973. 

2D.  J.  Pastine,  R.  R.  Bernecker,  and  R.  J.  Bauer,  "Theoretical  Relationship 
between  Initiating  Shock  Pressure  and  Porosity  in  Secondary  Explosives," 

Fourth  International  Conference  on  High  Pressure  (AIRAPT),  25-29  Nov  197I+,  in 
Kyoto,  Japan,  to  be  published  in  the  proceedings. 
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E^  calculated  from  the  critical  initiation  pressures,  Pc>  The  listed  porosity, 
p,  was  calculated  as  p * 1 - P/Pr^*  where  p is  the  initial  explosive  density 
and  is  the  theoretical  maximum  (voidless)  density  (TMD)  of  the  explosive. 
The  significance  of  the  results  in  the  table  are  discussed  below. 

It  is  emphasized  that  the  porosity/sensitivity  concept  requires  at  least 
one  experimental  sensitivity  datum  point  and  one  nonreactive  Hugoniot  equation 
of  state  for  the  explosive  before  it  can  be  applied. 

II.  OBJECTIVES  OF  THIS  PROGRAM 
The  purposes  of  this  program  of  analysis  were  three-fold: 

(a)  to  determine  the  validity  of  the  recently  developed  porosity 
concept2  derived  for  TATB-like  explosives  when  used  in  the  SSGT, 

(b)  to  determine  whether  this  concept  could  be  extended  to  other 
explosives,  and 

(c)  to  predict  the  sensitivity  of  a given  explosive  at  any  density 
from  60%  of  the  theoretical  maximum  density  up  to  the  TMD  with  a minimum  of 
experimental  data. 

The  explosives  discussed  in  this  report  are:  PETN,  TNT,  DATB,  TATB, 

DATB/Zytel  95/5 > and  DATB/Zytel  90/10.  The  required  Hugoniot  equations  of 
state  were  obtained  as  follows:  for  PETN  from  reference  3;  for  TNT,  DATB,  and 

TATB  from  reference  U;  the  values  for  DATB/Zytel  mixtures  were  assumed  to  be 
the  same  as  for  DATB  alone.  The  SSGT  values  were  obtained  from  reference  1. 
Using  equations  A-l,  A-5,  and  B-l,  E^  was  calculated  for  each  explosive  at 
each  porosity.  The  calculated  critical  pressures,  particle  velocities,  and 
thermal  energies  required  for  initiation  of  these  explosives  are  given  in 
Table  1 along  with  the  loading  pressures  and  calculated  densities  and  porosities 
of  the  explosives.  The  calculated  values  of  E^  for  each  explosive  were  examined 

and  it  was  found  that  is  in  fact  approximately  a constant  for  each  explosive 

for  all  observed  porosities.  E^  was  then  averaged  for  each  explosive  and, 

using  the  relationships  given  in  Appendix  B,  the  critical  initiation  pressure 
as  a function  of  porosity  was  calculated  and  plotted  for  each  explosive.  These 
plots  are  shown  in  Figures  1 through  6. 

3 J.  Roth,  "Shock  Sensitivity  and  Shock  Hugoniots  of  High  Density  Granular 
Explosives,"  Fifth  Symposium  on  Detonation,  18-21  Aug  1970. 

*N.  L.  Coleburn  and  T.  P.  Liddiard,  Jr.,  "Hugoniot  Equations  of  State  of  Several 
Unreacted  Explosives,"  J.  Chera.  Phys.  , UU_,  5,  1 Mar  1966. 
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III.  DISCUSSION  AND  CONCLUSIONS 

The  theoretical  predictions  of  the  critical  pressures  were  in  good  agreement 
with  experimental  results  with  certain  exceptions.  These  exceptions  are 
discussed  below.  The  agreement  of  theoretical  and  experimental  values  for  the 
plots  is  within  10$.  For  some  explosives,  the  agreement  is  much  better  than 
this. 

There  are  a number  of  possible  sources  of  error,  any  one  of  which  could 
be  responsible  for  the  disagreement  between  the  theoretical  and  experimental 
values  for  the  critical  pressure: 

(a)  Errors  in  determining  the  density  of  the  explosive  would  have  a 
significant  effect  on  the  critical  pressure,  especially  at  densities  near  the 
TMD. 

(b)  Batch-to-batch  variation  in  particle  size  and  chemical  purity 
could  cause  batch-to-batch  variation  in  sensitivity. 

(c)  The  nonreactive  Hugoniots  were  measured  at  low  values  of  particle 
velocity  and  linearized.  The  true  nonreactive  Hugoniot  is,  however,  nonlinear. 
Unfortunately,  the  sparsity  of  relevant  data  makes  it  difficult  to  determine 

| the  amount  of  nonlinearity.  It  would,  however,  have  some  effect  on  the  critical 

| pressure. 

I 

jj  (d)  The  porosity  concept  assumes  a uniform  distribution  of  voids. 

3 As  the  porosity  goes  to  zero  in  certain  explosives,  this  assumption  becomes 

\ less  and  less  valid.  The  nonuniform  distribution  of  voids  causes  the  sensitivity 

| to  decrease;  that  is,  the  actual  critical  pressure  is  significantly  greater 

■ than  the  predicted  value.  This  problem  is  demonstrated  on  the  TNT  curve,  for 

j,  porosities  less  than  0.017. 

For  TATB  at  porosities  above  0.075,  "the  predictions  are  good,  but  at  lower 
values  of  porosities  (higher  densities),  the  curve  differs  significantly  from 
; experimental  values.  Since  the  porosity  concept  had  been  developed  for  TATB  and 

\ the  only  material  which  did  not  have  good  predictability  in  the  SSGT  was  TATB, 

; an  apparent  contradiction  existed. 

| The  available  SSGT  calibration  curve  relating  pressure  to  barrier  thickness 

} is  linear  for  barrier  thicknesses  greater  than  2.5  mm.  It  is  non-linear  and  not 

well  defined  for  lesser  thicknesses.  Unfortunately,  the  barrier  thicknesses  for 
j the  50#  functioning  points  of  TATB  at  the  densities  of  1.8k0  and  1.887  (Table  l) 

were  smaller  than  2.5  mm.  Thus,  the  estimate  uf  the  critical  pressure  at  these 
two  densities  was  subject  to  considerable  error  and  is  most  likely  the  cause  of 
the  poor  agreement  found  at  these  densities. 


^ 5D.  Price  and  T.  P.  Liddiard,  .Jr.,  "The  Small  Scale  Gap  Test:  Calibration  and 

s Comparison  with  the  Large  Scale  Gap  Test,"  NOLTR  66-87,  7 Jul  I966. 
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On  the  basis  of  the  good  agreement  between  the  theoretical  prediction  of 
the  critical  pressure  and  the  experimental  data*  and  the  fact  that  the  few 
differences  can  be  reasonably  explained,  these  studies  support  the  validity  of 
the  recently  developed  porosity  concept  for  most  secondary  explosives.  It  is 
also  apparent  from  the  fit  of  the  critical  pressure  versus  porosity  curve  of 
Figures  5 and  6 that  the  concept  holds  true  for  explosive  compositions  desensitized 
with  nonreactive  materials,  e.g.,  binders,  lubricants,  pelletizing  agents. 

Previously,  to  estimate  the  SSGT  shock  sensitivity  cf  an  explosive  material 
at  a density  other  than  at  a testing  value,  one  would  interpolate  or  extrapolate 
the  data  from  the  two  closest  points  or  from  a curve  fitted  to  all  of  the 
available  data.  These  methods  are  not  based  on  the  operative  hydrodynamic  and 
thermodynamic  parameters,  and  are  therefore  of  questionable  accuracy. 

The  porosity/sensitivity  concept,  on  the  other  hand,  enables  U3  to  generate 
other  Hugoniot  equations  of  state  and  to  predict  sensitivities  using  a single 
Hugoniot  and  a single  sensitivity  point.  These  latter  values  are  often  available 
for  many  explosives.  This  method  has  a distinct  advantage  because  no  further 
interpolation  or  extrapolation  of  the  data  is  necessary,  and  one  has  a much 
greater  degree  of  confidence  in  the  result  than  can  be  obtained  from  the  SSGT 
data  alone. 

To  make  estimates  of  explosive  system  reliability  and/or  safety,  one  must 
know  the  critical  pressure  required  to  initiate  the  explosive.  Since  the 
material  we  are  to  initiate  is  the  explosive,  it  is  more  useful  to  know  the 
critical  pressure  in  the  explosive  than  the  critical  pressure  in  the  PMMA 
attenuator  of  the  SSGT.  Utilization  of  the  porosity/sensitivity  concept  permits 
the  computation  of  the  desired  parameter — the  critical  pressure  in  the  explosive  — 
if  one  Hugoniot  equation  of  state  at  a given  porosity  is  available. 

The  porosity  concept  appears  applicable  to  desensitized  explosive 
compositions.  It  would  be  valuable  to  have  a theoretical  relationship  between 
the  percentage  of  diluent  in  the  composition  and  the  critical  thermal  energy. 

One  possible  approach  would  be  to  treat  the  diluent  as  "solid  holes".  This 
work  would  complement  the  current  concept. 
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TABLE  1 POROSITY /SENSITIVITY  DATA  (Continued) 


Explosive  Loading 
Pressure 
(kpsi) 


Density  Porosity 
(gm/cc) 


c 

(Kber ) 


Notes 


(mm/usec)  (Joules /min'3) 


1* 

1.231 

.3310 

23.56 

.8711 

.1+8217 

1* 

1.233 

• 3299 

23.78 

.871+0 

.1+8588 

1+ 

1.255 

.3179 

23.93 

.8573 

.1+7257 

8 

1.339 

.2723 

29.61 

.8889 

.52531+ 

8 

1.365 

.2582 

28.18 

.81+00 

.1+7227 

10 

1.1+55 

.2092 

31+.83 

.8653 

.50681+ 

13 

1.  UU2 

.2163 

35.26 

.881+2 

.52929 

16 

1,511+ 

.1772 

37.1+8 

.81+76 

.1+8327 

16 

1,518 

.1750 

36.26 

.8265 

.1+5886 

1 6 

1.655 

.1005 

5I+.I+6 

.9160 

.52811 

32 

1.662 

.0967 

51.56 

.8757 

.1+7638 

32 

1.665 

.0951 

51.23 

.8687 

.I+669I+ 

32 

1.676 

.0891 

51+.00 

.8881 

.1+81+18 

32 

1.738 

.0551+ 

63.39 

.9179 

.1+8125 

37 

1.701 

.0755 

56.61+ 

.8896 

.1+7217 

50 

1.732 

.0587 

6I+.1I+ 

.9322 

.50311+ 

6k 

1.763 

.01+18 

76.69 

1.0155 

.58709 

6k 

1.775 

.0353 

71+.30 

.9790 

.52708 

Average 

.1+9151 

k 

1.519 

.2130 

1+3.19 

.91+27 

.67332 

8 

1.61+5 

.11+77 

60.06 

1.0061 

.77597 

8 

1.727 

.1052 

83.29 

1.1276 

.97953 

16 

1.762 

.0870 

93.21 

I.16U7 

1.01+023 

32 

1.81+0 

.01+66 

: ’+2.1+7 

1.3969 

1.55119 

6k 

1.887 

0223 

312.80 

C.  161+3 

1+.1+0771+ 

Average 

.86721+ 

NOTES:  * Items  not  used  in  computing  average 
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TABLE  1 POROSITY/SENSITIVITY  DATA  (Continued) 


Explosive 


DATB/ZYTEL 

95/5 


DATB/ZYTEL 

90/10 


Loading 

Pressure 

(kpsi) 


Density 

(gm/cc) 


Porosity 


c 

(Kbar ) 


(mm/Usec)  ( Joules /mm^) 


1* 

1.192 

.3522 

28.26 

1.001+1 

.62939 

1* 

1.210 

.31+21+ 

28.79 

.9969 

,62656 

8 

1.358 

.2620 

37.01* 

.991*5 

.66317 

8 

1.366 

.2576 

37.02 

.9862 

.65321 

16 

1.529 

.1690 

1+8.86 

.981+5 

.65507 

16 

1.534 

.1663 

1+6.60 

.9501+ 

.60806 

32 

1.657 

.0995 

59.1+2 

.9681 

.591+1+1 

32 

1.661 

.0973 

65.32 

1.0258 

.67290 

Avei age 

.63785 

1* 

1.167 

.3658 

30.23 

1.0687 

.701+19 

8 

1.31+2 

.2707 

38.76 

1.0380 

.72066 

16 

1.512 

.1783 

52.53 

1.01+93 

.75063 

32 

1,617 

.1212 

61+ . 60 

I.0691 

.75727 

6k 

1.676 

.0691 

77.69 

1.1*07 

.82780 

Average  .75217 


^ar.  jt&i-  ’ -i&i  r&L  *iiiiWJ;*  ji  . ,i*5f-  fiiliA  -.j,  i = ^ -ii.wt.iJ  jJft..  .ijj..  j.j  jlliC  2 *2j)lifL^  k Is.  1H1J  J.i-:>  .^lidUai.it*!#  rt^ya^iJ  ;:*i.  ..t*  - ..  £^.i  ],(JjiJJji.j.^>  utf 
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Et=  0.06022  JOULES/MM3 

X'S  ARE  EXPERIMENTAL  POINTS 

SOLID  LINE  CURVES  ARE 
CALCULATED  VALUES 


FIG.  1 SENSITIVITY  VS  POROSITY  CURVE  FOR  PETN 
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POROSITY 

FIG.  4 SENSITIVITY  VS  POROSITY  CURVE  FOR  TATB 
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APPENDIX  A 

CALCULATION  OF  SHOCK  HUGONIOTS 

1.  The  nonreactive  shock  Hugoniot  for  seme  given  secondary  explosive  at  a 
particular  density  is  calculated  according  to  the  porosity  oonceptA“l  in  the 
following  way: 


up2(P,p)  * up2<p,0)  * ^ 


TMD 


1 - P 


(1) 


where  u^  is  the  particle  velocity  of  the  explosive, 

P is  the  pressure  in  the  explosive, 
p is  the  porosity  of  the  explosive,  and 

PTMD  is  the  theoretical  maximum  density  (TMD)  of  the  explosive. 


U = C + Bu 
s p 


and 


P = p u U , 
0 o s 


(2) 


(3) 


where  U is  the  shock  velocity  in  the  explosive,  C and  B are  the  constants  of 
s 

the  linearization  equation  (C  being  the  sound  speed  intercept  and  B,  the  slope) 
and  pQ  is  the  density  at  which  the  nonreactive  shock  Hugoniot  was  determined. 

p 

2.  Inserting  equation  (2)  into  (3)  and  solving  for  u , we  obtain, 

P 


1 + 2 


PB 

po°E 


■Vi* 


PB 

PoC' 


(U) 


where  p^  is  the  porosity  associated  with  p , by 


Po  ” 1 “ Po/pTMD  ' 

2 

Inserting  equation  (1+)  into  (l),  substituting  pQ  for  p,  and  solving  for  (P,0), 
we  get 


A-l 
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between  Initiating  Shock  Pressure  and  Porosity  in  Secondary  Explosives," 
published  in  the  Proceedings  of  the  Fourth  International  Conference  on  High 
Pressure  (AIRAPT),  Kyoto,  Japan,  Mar  1975. 
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PB 

poC' 


P, 


2 Pi 


TMD 


1 - P, 


(5) 


With  equation  (5)»  equation  (i)  may  be  solved  for  any  porosity. 

A 2 

3.  We  now  make  use  of  shock  impedance  matching  theory  ~ to  solve  for  the 
critical  pressure  in  DATB  at  a density  of  1.775  gm/cc:  PTMD  = 1.8U  gm/cc, A-3' 

B ■ l.892,A_1*  C = 2.UU9  mm/ysec,A_!+  Pq  = 1.780  gm/cc ,A~^  pQ  = 0.0326,  and  p = 0.353 

Inserting  these  numbers  into  equations  (l)  and  (5)  and  plotting  P as  a function 
of  up,  we  set  the  solid  line  in  Figure  A-l. 


The  Hugoniot  equations  for  PMMA  are 


A- 5 


P = pu  U 
P s 


U = 
s 


2.7228  + 1*.o667u  - 10.9051u  2 + 10.6912u  3, 

P P P 

for  0.03  mm/ysec<u  <0.5363  mm/ysec 
P 

2.56l  + 1.595^,  for  Up>0.5363  mm/ysec 


where  P is  the  pressure  in  the  PMMA, 


(6) 

(7) 

(8) 


Ug  is  the  shock  velocity  in  the  PMMA, 

^ is  the  particle  velocity  in  the  PMMA,  and 
p is  the  original  density  of  the  PMMA  (1.185  gm/cc). 

Plotting  P as  a function  of  u^  we  get  the  dashed  line  in  Figure  A-l. 

At  a density  of  1.775  gm/cc,  the  SSGT  sensitivity  of  DATB  is  8.882  DBg,A-3 
which  corresponds  to  a shock  pressure  of  60.22  Kbar  in  the  PMMA.A“6  A PMMA 
Hugoniot  reflected  about  the  vertical  line  passing  through  60.22  Kbar  on  the  PMMA 
Hugoniot  is  shown  in  Figure  A-l. 


A -2, 


A- 3 


A-k 


A- 5 


A- 6. 


G.  E.  Duvall  and  G.  R.  Fowles,  "Shock  Waves,"  High  Pressure  Physics  and 
Chemistry,  2_,  Chapter  9,  Academic  Press,  1963. 

o'.  N.  .Ayres,  L.  J.  Montesi,  and  R.  J.  Bauer,  "Small  Scale  Gap  Test  (SSGT)  Data 
Compilation:  1959-1972  Volume  I Unclassified  Explosives,"  N0LTR  73-132, 

26  Oct  1973. 

N.  L.  Coleburn  and  T.  P.  Liddiard,  Ti . , "Hugoniot  Equations  of  State  of  Several 
Unreacted  Explosives,"  J.  Chem.  Pnys.,  No.  5,  1 Mar  1966. 

J.  0.  Erkman,  D.  J.  Edwards,  A.  R.  Clairmont,  Jr.,  and  D.  Price,  "Calibration 
of  the  N0L  Large  Scale  Gap  Test;  Hugoniot  Data  for  Polymethyl  Methacrylate," 
N0LTR  73-15,  h Apr  1973. 

D.  Price  and  T.  P.  Liddiard,  Jr.,  "The  Small  Scale  Gap  Test:  Calibration  and 
Comparison  with  the  Large  Scale  Gap  Test,"  N0LTR  66-87,  7 Jul  1966. 

A-2 


N-JWC/WOL/TR  75-67 


The  reflected  PMMA  Hugoniot  intersects  the  nonreactive  shock  Hugoniot  for 
DATE  at  1.775  gm/oc  at  a pressure  of  7^.30  Kbar  and  a u^  of  0.979  mm/ysec. 

These  are  respectively  the  critical  pressure  and  critical  particle  velocity  of 
DATB  at  an  original  density  of  1.775  gm/cc.  This  result  was  obtained  by  a 
graphical  solution  in  Figure  A-l,  but  was  performed  by  a computer  program  as  an 
iterative  solution. 
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APPENDIX  B 

CALCULATION  OF  CRITICAL  ENERGY  DENSITIES 

The  specific  critical  thermal  energy  per  unit  volume  is  calculated 
with  the  relation, 


PTMDC 


JM. 


1"n  I 2(1  - nB)2  2(1  - nB)‘ 


-*[ 


(1  - nB)-1  + ln(l  - TIB)  - 1 


(1) 


where  E.  is  the  critical  thermal  energy  per  unit  volume, 

IT  is  a function  of  porosity:  II  = p/(l  - p), 

n is  a function  of  the  specific  volume  behind  the  initiating  shock  wave: 

n ■ 1 - y/y±t 

vi  " 1//pTMD’ 

p is  the  initial  porosity  (the  fraction  of  the  initial  volume  which 
consists  of  holes,  and 

B and  C are  as  defined  in  the  text. 

For  the  example  of  the  DATB  at  a density  of  1.775  gm/cc , E^  is  0. 52708  Joules/mm  . 

To  generate  the  critical  pressure  versus  porosity  plots,  ve  made  use  of  tho 
following  equation: 


n 2d  - nB)g  Vi(1  - 111 


?[ 


2(1  - nB)‘ 


(1  - nB)-1  + ln(l  - nB) 


■*] 


(2) 


where  is  the  average  critical  thermal  energy  per  unit  volume  for  a given 
explosive  (except  for  the  asterisked  items  in  the  "Notes"  column  of  Table  l). 

Equation  (2)  relates  the  porosity  to  the  specific  volume  behind  the 
initiating  shock  wave.  In  order  to  relate  the  porosity  to  the  critical  pressure, 
the  latter  volume  was  used  to  find  the  associated  shock  pressure  by  means  of  the 
theoretically  determined  shock  Hugoniot  for  the  porosity  under  consideration. 
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